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Abstract 

Elastic collisions in the transplanckian region, where the center-of-mass energy is 
much larger than the fundamental gravity mass scale, can be described by linearized 
general relativity and known quantum-mechanical effects as long as the momentum 
transfer of the process is sufficiently small. For larger momentum transfer, non-linear 
gravitational effects become important and, although a computation is lacking, black- 
hole formation is expected to dominate the dynamics. We discuss how elastic trans¬ 
planckian collisions can be used at high-energy colliders to study, in a quantitative 
and model-independent way, theories in which gravity propagates in flat extra dimen¬ 
sions. At LHC energies, however, incalculable quantum-gravity contributions may 
significantly affect the experimental signal. 
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1 Introduction 


The hypothesis that the characteristic energy scale of quantum gravity lies just beyond the 
electroweak scale hnds its motivation in the hierarchy problem and its operative realization in 
the existence of extra spatial dimensions 1^,0. If this hypothesis is verihed in nature, future 
experiments at the LHC have the ground-breaking opportunity of testing the dynamics 
of gravitation in the high-energy and quantum regimes. Unfortunately, since quantum- 
gravity dynamics are still unknown, systematic theoretical predictions cannot be made. It 
is then important to identify certain collider topologies and kinematical regions that allow a 
description based only upon general principles and not on model peculiarities, and which can 
be used in the experimental program as a handle to demonstrate the gravitational nature of 
the interaction. Single graviton emissions satisfy these requirements and allow for a model- 
independent study of the propagation of gravitational forces in extra spatial dimensions p, 

On the other hand, model-dependent experimental signals, as those originating from contact 
interactions p, P or, more likely, the observation of a specihc spectrum of string Regge 
excitations [^] will be very important to discriminate among various models and to guide 
theoretical research. 

The model-independent experimental signals discussed in refs. P P are based upon 
interactions at energies below the quantum gravity scale Md (the analogue of the Planck 
mass for a D-dimensional theory), and their theoretical description fails as we approach 
Mq. However, there is another kinematical regime that can be tackled in a fairly model- 
independent fashion: the transplanckian region -y/s S> Md- This is the subject of this paper. 
A pecularity of hadron machines is that, because of the composite structure of the proton, 
one can study, in the same environment, parton collisions at different center-of-mass energies. 
Consequently, if ~ TeV experiments at the LHC could probe the cisplanckian (graviton 
emission, contact interactions), the planckian (the core of the quantum-gravity dynamics), 
and the transplanckian regions. 

To understand the nature of transplanckian collisions, we hrst identify the relevant scales 
(see ref. P). Let us consider general relativity in D = 4 -|- n dimensions with (generalized) 
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Newton’s constant Gd- For n flat dimensions compactified in a volnme 14i, Gd is related 
to the usual Newton’s constant by Gd = VnGN- To better elucidate the relation between 
the transplanckian and the classical limit, we reinstate for a moment the correct powers of 
h and c. The relation between Gd and the D-dimensional Planck scale is 




Gn = 


{2Tr)^-^K 
4c«-iM]^+2 ’ 


( 1 ) 


where the proportionality constant follows from the conventions of ref. |0]. Notice that Gd 
has dimensions [Gd] = {i = length, E = energy, t = time). Starting from Gd, 

we can construct the Planck length 


GDh\^ 

) 


( 2 ) 


This is the distance below which quantum gravity effects become important. 


Using the center-of-mass energy of the collision, we can construct the length 

AttHc 


Xb — 


V~s ■ 


( 3 ) 


This is the de Broglie wavelength of the colliding particles, which characterizes their quantum 
length scale. 


Combining Gd and a/s, we can form the Schwarzschild radius of a system with center- 
of-mass energy a/s [|| 


Rs 



8r(:f!) 

(n + 2) 


1 

n+1 



( 4 ) 


This is the length at which curvature effects become significant. In the limit h —>■ 0, with 
Gd and ^/s hxed, Md vanishes, showing that classical physics correspond to transplanckian 
(macroscopically large) energies. Moreover, in the same limit, Rs remains finite, while the 
two length scales Ap and Ap go to zero. Therefore, the transplanckian regime corresponds 
to a classical limit in which the length scale Rs characterizes the dynamics. 


a/s S> Md Rs 3> Ap S> Ap. 


( 5 ) 
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For instance, by simple dimensional analysis and analiticity in Gd, we expect that the 
scattering angle for a collision with impact parameter b is 6 = {Rs/br+^. 

This behavior shows that by increasing the energy we can obtain a hnite angle scattering by 
going to larger b, and therefore further suppressing short distance quantum gravity effects. 

This property of gravity should be contrasted to what happens in an ordinary gauge 

theory, where spin-1 particles mediate the force. In (4 -|- n)-dimensional space-time the 

gauge coupling has dimension [e^] = The analogues of \p and Rs are given 

respectively by M” = , A” = /{Tic) and R/'^^ = In the h —0 limit, with 

and ^/s hxed, we hnd that \e ^ oo (we focus for simplicity on n > 0) indicating that 

quantum fluctutations are not suppressed at any hnite length. Indeed the same dimensional 
argument we used before gives a classical scattering angle 6 ~ which remains 

hnite as y/s oo provided 6 —0. But b < Xg is the regime where quantum ehects are 
unsuppressed, so we conclude that there is no classical limit. The diherent properties of spin- 
2 and spin-1 exchange that are elucidated by this simple dimensional analysis are the same 
that render the eikonal approximation consistent for the hrst and inconsistent for the second. 
Basically it is because energy itself plays the role of charge in gravity: when the energy is 
large, gravity gives sizeable ehects also at large distance where quantum huctuations of the 
geometry are irrelevant. 

Having established that scattering at transplanckian energies is described by classical 
physics, one also realizes that the corresponding amplitude should only be calculable by 
a non-perturbative approach. This is always the case for semiclassical amplitudes (like 
instantons for example), for which the classical action S is much larger than h, and we 
cannot perturbatively expand exp(—S'/h). In order to tackle the non-perturbative problem, 
we will further restrict the range of distance scales under investigation, and only consider 
collisions with impact parameter b much larger than Rs- In this regime, the curvature is 
small, the metric is nearly hat, and we can work in the limit of weak gravitational held, 
neglecting non-linear ehects of the graviton couplings. The non-perturbative interactions 
between the high-energy colliding partons and the weak gravitational held can be computed 
using the eikonal approximation, which can be trusted at small scattering angle. From a 


4 



field-theoretical point of view, this approximation corresponds [Q to a resummation of an 
inhnite set of Feynman diagrams which, at each order in perturbation theory, give the leading 
contributions to the forward scattering. The eikonal approximation has been employed for 
4-dimensional gravity and for string theory [[^ (see also ref. []K|). The results for the 
scattering amplitude agree with those obtained []^ by solving the Klein-Gordon equation 
for a particle propagating in the classical gravitational shock wave produced by the other 


particle |I^, or by reformulating quantum gravity as a topological held theory [|I^. All these 


different approaches give equivalent results 10 


In sect. § we will derive and discuss the scattering amplitude in the eikonal approximation 


for 4-dimensional colliding particles (living on a 3-brane) and D-dimensional gravity (see 
also ref. [T^). As explained above, the process is essentially classical, like the motion of two 
(relativistic) stars coming within a close distance. However, an important difference is that 
the gravitational potential in the non-relativistic analogue of our process decreases with the 
distance r as H(r) oc As discussed in sect. | this will have signihcant consequences, 

and it will lead to quite different results from the case of 4-dimensional Newtonian potentials. 
A novel feature is the emergence of a new length scale 


be. oc 


V hc^ ) 


( 6 ) 


Notice that be cannot be dehned in 4 dimensions (n = 0). In the limit h —>■ 0, we hnd that 
be goes to inhnity and therefore the classical region extends only up to length scales of order 
be- For impact parameters larger than be, the process is no longer determined by classical 
properties. However, since quantum gravity effects are always negligible, ordinary relativistic 
quantum mechanics is sufficient to compute scattering amplitudes within our approximation. 


As previously discussed, we are limited by our computational ability to work at 6 S> Rs- 
As the impact parameter approaches the Schwarzschild radius (and our expressions become 
unreliable), we expect to enter a (classical) regime in which the gravitational held is strong, 
non-linear effects are important and, plausibly, black holes are formed. This regime, discussed 
in refs. |jT^ and (see also ref. |^), is very exciting from the experimental point of view. 


but more arduous from a theoretical point of view. While in the eikonal region we are 
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able to make quantitatively reliable predictions with controllable expansion parameters, the 
black-hole production cross section can only be estimated with dimensional arguments and 
more precise information should wait for numerical simulations in D-dimensional general 
relativity. 


Black-hole production results in a black disk of radius Rs in the cross section, since any 
initial state with b < Rs is absorbed by the black hole (this point of view has been criticized 
in ref. |^^). Collisions with an impact parameter b smaller than the Schwarzschild radius 
do not directly probe physics at the energy scale 1/6. Ultraviolet dynamics and quantum 
effects are screened by the black hole. 


The physical picture we have been describing holds if there is a clear separation of the 
various length scales and if no other new dynamics intervene at distances between Rs and 
be- As we will discuss in this paper, this is not necessarily the case for particle collisions 
at the LHC. Given the present lower limits on M^, the center-of-mass energy at the LHC 
can be only marginally in the transplanckian region, and a clear separation between the 
quantum-gravity scale \p and the classical scale Rs is not completely achieved, endangering 
our approximations. The problem is particularly apparent in string theory, which contains 
a new scale, the string length 

As = (7) 


where a' is the string tension. In weakly-coupled string theory. As is larger than Ap and 
string effects can set in and modify gravity in an intermediate regime between the eikonal 
elastic scattering and black-hole formation |^. These are interesting effects, as the LHC can 
directly probe the nature of string theory [^]. However, model-independent calculability 
of the elastic channel is lost. In this paper we will consider only the pure transplanckian 
gravitational scattering, assuming that quantum-gravity effects are small. 


The relevance of our considerations for experiments at the LHC and other future colliders 
is discussed in sect. 3. We will compute the elastic cross section in the transplanckian eikonal 
regime, 

Afp»/\/s-Cl, — f/s<Cl, (8) 
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where t is the Lorentz-invariant Mandelstam variable, and —t/s is a measure of the scattering 
angle in the center-of-mass frame. The experimental signal will be discussed and compared 
with the expected Standard Model (SM) rate and with black-hole production. 


2 Scattering in the Eikonal Approximation 

2.1 The Eikonal Amplitude 


We are interested in the elastic scattering of two massless 4-dimensional particles (living on 
a 3-brane) due to a D-dimensional weak gravitational held in the kinematic limit of eq. (^. 
We will focus on the case in which the two colliding particles are different, but the case of 
identical particles is analogous. In the eikonal approximation, we consider the inhnite set of 
ladder and cross-ladder Feynman diagrams which give the leading contribution to forward 
scattering. 


Since we are working in the approximation of small momentum transfer, we can compute 
the coupling between the colliding particle $ and the small huctuation around the hat metric 
(the graviton), taking $ to be on-shell: 

(9) 

Here is the energy-momentum tensor and is the 4-momentum of the particle <h. Notice 
that eq. (^ is equally valid for colliding particles of any spin. 

In the <F propagators we keep only the leading terms in the momentum transfer q: 


i 

{p + qY 


i 

j/+ (/ 


- for bosons, 

2pq 


i|4>(p))(4>(p)| 

2pq 


for fermions. 


( 10 ) 


Since the factor ]/ = |4>(p))(<F(p)| in eq. m is used to reconstruct the on-shell vertices of 
eq. (^) in the ladder diagrams, we observe that not only the interaction vertices but also 
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Figure 1: The dominant Feynman diagrams contributing to elastic scattering of $ parti¬ 
cles in the eikonal approximation at tree (a) and one-loop (6) level. Wavy lines represent 
the exchange of D-dimensional gravitons, and the dashed line represents the cut from the 
physical-region singularity. 


the propagators of bosons and fermions are identical in the eikonal approximation. This 
universality of bosons and fermions is not surprising since, as discussed in the introduction, 
we are performing a semi-classical calculation, and the spin of the <F particle should not 
matter. 


The tree-level exchange of a D-dimensional graviton in the diagram of £g. gives the 
scattering amplitude 


A 


Born I, 


-i) = 




d^OT , s ( —t 

=7r2F(l-u/2) 


r-l 


t-Qr 


Ml, 


Ml, 


( 11 ) 


where qt is the momentum transfer in the extra dimensions. Although eq. ( 0 ) is a tree- 
level amplitude, it is divergent due to the inhnite number of (extra-dimensional) momentum 
conhgurations of the exchanged gravitons, allowed by the non-conservation of momenta 
transverse to the brane. Divergences have been subtracted using dimensional regularization, 
considering non-integer n, but the regularization prescription is not important since the 
eikonalization will consistently eliminate any ultraviolet sensitivity. Basically this is because 
the eikonalization selects in eq. (pd]) only the partial waves with large angular momentum. 
These do not depend on the local counter-terms as they are determined by the hnite (calcu- 











lable) non-analytic terms in —t, which depend only on infrared physics |^. These terms are 
proportional to (—ln(—t) for even n, and to a/^(— for odd n. 


Let us now turn to the one-loop amplitude. The leading term in the limit of small 
momentum transfer is given by a Feynman diagram with a cut arising from a physical-region 
singularity in which the intermediate $ particles are on their mass shells. There is only one 
such diagram at one loop, and it is drawn in £g. B). All other diagrams are subleading. 
For instance, the one-loop t-channel exchange of two gravitons G coupled via quartic 
vertices is suppressed^ by a power of t/s. Diagrams with one vertex and two 

vertices simply vanish (as is evident by working in the de Bonder gauge). 


The diagram of fig. |J) with on-shell intermediate states can be calculated using the 
Cutkosky rule. Its imaginary part is 1/2 times the discontinuity across the cut, which is 
obtained by replacing the $ propagators in eq. m with 27ri6{2pq). If pi and p 2 are the 
momenta of the incoming particles, we obtain 


•^1—loop( Q ) 


d?ki 

W) 




d^k 


rkr^k'^ 


-k^{q- kf - k'S 


r'ABorn(^_L)“^Born[(?± k±} 


{2'Kif‘ 5{2pik)5{2p2k) 


( 12 ) 


The delta functions from the Cutkosky rule had the effect of reducing the 4-dimensional 
integral into a 2-dimensional integral over momenta perpendicular to the beam (and along 
the brane). Notice that, in the eikonal limit, the momentum transfer q is mainly in the 
perpendicular direction and therefore, we can use the approximation 


j. 2 2 

t = q ^ -q±. 


(13) 


Equation (|T^ is merely a convolution of the Born amplitude. It is then convenient to 

perform a Fourier transform of the amplitudes with respect to the transverse momentum qj_. 

This amounts to trading qj_ with its conjugate variable bj_, the impact parameter. In impact 

^Actually, in the spirit of the eikonal resummation, this diagram should be compared to the tree-level 
graviton exchange, but an analogous suppression factor would nevertheless appear. 
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parameter space eq. ([T^ ) becomes a simple product. Summing eq. ( p!TD and eq. and 

transforming back to momentum space we can recast the loop expansion in the form 


^Born(gi) + A-ioop(<?l) + • • • = -2is J i^ix 


2 % +••• 




(27r 


(14) 

(15) 


The combinatorics of higher-order loop terms is such that one can resum all terms in eq. ([T^) 
to obtain 

( 16 ) 


Aik = - 2 m [ (e‘» - l) 


We have recovered the known result @ that the eikonal scattering phase Xi function of the 
impact parameter plane coordinates b±, is the 2-dimensional Fourier transform of the Born 
amplitude in the direction perpendicular to the beam. 

One crucial feature of the eikonal procedure is that x(^) for 6 7 ^ 0 depends only on the 
calculable, ultraviolet hnite, terms in A-Bom- The ultraviolet divergent terms correspond to 
deZfa-function singularities localized at 6 = 0. While these terms are obviously irrelevant 
for scattering processes where the initial states are prepared in such a way that they have 
no overlap at & = 0 , we will also dehne our full amplitude in eq. m by neglecting these 
terms. We believe this to be a consistent procedure. Indeed, it would not make any sense 
to include ^-functions in the exponent of eq. (pT]). From a physical point of view, we expect 
these singularities to be softened by the fundamental theory of gravity, at some hnite, but 
small, impact parameter 6 ~ Ap (or, more likely, b ~ As). So the presence of these terms 
is just evidence that there should be corrections to our simple eikonal expression coming 
from the underlying theory of quantum gravity. We will discuss these effects in sect. ^]5. 
For the moment we limit ourselves to note that these short-distance contributions should 
plausibly give rise to C>(Ap) corrections to the cross section, while, as we will see shortly, the 
long-distance eikonal amplitude gives a cross section that grows with a power of y/s, thus 
dominating at large energies. 
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We can now calculate explicitly the eikonal phase y in eq. ( 0 - One hnds (see appendix) 

(17) 


ttS W(l-n/2)s 

X = - 777 ^ 1+2 - I dqq Mqb) = ( - 




br = 


/o 


(dvr) 2 ^sr(n/2) 


2 M ])+2 


l/n 


(18) 


Here g = |g_i_| ~ \/^, b = | 6 _l| and Jo is the zero-th order Bessel function. Notice that by 
inserting this result for y in the integral in eq. ( 0 ) we obtain an ultraviolet hnite result. This 
is so, even though the contributions to eq. (P!B| ) from the individual terms in the expansion 
of = 1 + iy + ... are ultraviolet divergent, corresponding to the fact that each individual 
Feynman diagram in the ladder expansion is ultraviolet divergent but the complete sum is 
hnite. Moreover, since y oc 6 “"', the integrand in eq. (|I 6 D oscillates very rapidly as 6 —> 0 , 
showing that the ultraviolet region gives but a small contribution to the amplitude. 

Replacing eq. (0) into eq. (|I 6 D , we obtain the hnal expression for the eikonal ampli¬ 
tude |T7| 

Aeik = i7isblFn{bcq), (19) 


- 1 


( 20 ) 


^n( 2 /) = -«^ dxxJo{xy)(e 

where the integration variable is x = b/bc. The functions Fh, which can be written in terms 
of Meijer’s G-functions, are shown in £g. 0. 


At small y, the functions Fn{y) can be expanded as (see appendix) 

F2{y) 


Fsiy) 

F,{y) 

Fn{y) 


= 1 h -TT 

-ln_ + ,- 


= I 


j/^0 -T - 




-y 
,2 
4 


i—^e 4 -f — In?/ 


= i 


y^O -T 1- 


n 


- -y^T 1 - - 


16 


n 


for 77, > 4. 


( 21 ) 

( 22 ) 

(23) 

(24) 


These functions develop non-analytic terms in y'^ of the form y'^~‘^\ny for even n, and y'^~‘^ 
for odd n, which reproduce the non-analytic structure of M-Bom previously mentioned. Notice 
that Fh( 0 ) is hnite for n > 2. 


11 









0 2 4 


80246802468 


Figure 2: Plots of the function Fn{y) versus y for n between 2 and 6. The dashed line is the 
real part, the dash-dotted line is the imaginary part, and the solid line is the absolute value 
of the function. The bottom right panel plots the relative error of \FQ{y)\ compared to the 
asymptotic expressions, see Eqs. and as y —> 0 (dashed line) and y ^ oo (solid 
line). 
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For 3> 1 (<? ^ &c the stationary-phase approximation applies, yielding 


1 n + 2 

= —iri'^+^y "+i 

Fn{y) ?/>i- / — exp 

Vn + I 


-,(u + 1) (P 


( 26 ) 


The function oscillates around its asymptotic value given in eq. before converging 
to it at large y, as illustrated in the bottom-right panel of £g. |[ 


2.2 Cross Sections 


Let us now consider the differential cross section for elastic scattering. From eq. (0 we 
obtain 


We recall the dehnitions 

Ris = A„ 


dcTeik ,4 
—;— = 7r5_ 
dt 


71 + 2 
71 + 1 


Ml, 


Fn{bc\/^) 


bis = Br, 


71 + 2 


Mi. 


(26) 


(27) 


where the coefficients An and Bn are given in table 1. In terms of the geometric black-hole 


7 ri?|, eq. 

(p6D can be written as 

dcT^ik 

71 

Bn^^ aBH 


71 + 2 

W+1 

dt 

An ^ 

(n+2)2 



bI(Tbh / 


71(71+1) 


An S \ 

Ml) 



FnibcV^) 


(28) 


In the limit q ^ b^^, or equivalently —t/s <C the differential cross 

section becomes (for n > 2) 

71 

daeik [F(l - 2/n)]^Bn^^ (Jbh /,2 aStt 


dt 


4A 






[T{l-2/n)]‘^BlaBH f s W("+i) 


4:Ar. 


s \Ml 


(29) 


13 



























n = 

2 

3 

4 

5 

6 

An — 


^3 / 1 o \ 

2'"7r^“r(2M) 

2 

n+1 

0.83 

1.79 

3.01 

4.42 

5.96 


n-\-2 



Bn = 


2"-37r'^r (f 

)1 

2_ 

n 

0.50 

1.35 

2.51 

3.88 

5.41 


Table 1: The coefficients An and Bn entering eq. (^) and their numerical values for n 
between 2 and 6. 


In the opposite limit {q ^ we hnd 

do-eik _ (TsH / sxSrf 

dt {n + l)An s \ t) ' ^ ’ 

Notice how the two formulae valid in the two regimes match since, for —t ~ 6“^, eq. (|30|) 
parametrically agrees with eq. (^). Moreover, extrapolating eq. (^) beyond its range of 
validity towards small impact parameters (he. large angle, or t —> —s), we hnd a result that 
matches parametrically the black-hole production cross section derived from dimensional 
analysis. 


The integrated elastic cross can be derived directly from eq. (pTf): 


O'el 


J d%^_ (l + cos Rex) • 


(31) 


The total cross section, inclusive of elastic and inelastic channels, can be derived from the 
optical theorem 

o'tot = ^ ^ ~ cos Rex) • (32) 

In the absence of absorptive parts (Imx = 0) eq. (|3l|) and eq. (|3^ coincide. Black-hole 
production may be modelled by a large absorptive part (Imx 3> 1) at impact parameters 
smaller than Rs- 
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2.3 Physical Interpretation 


Before proceeding to the phenomenological analysis, we remark on the physical meaning of 
the expression of the eikonal amplitnde given by eq. (|^) and eq. (^). We can hrst consider 
a non-relativistic analogne of the problem nnder investigation by stndying the scattering 
of a particle with mass m, velocity n, and impact parameter 6 by a radial “gravitational” 
potential 


V{r) = 


GomM 


r-n+l 


(33) 


Here m and M are the non-relativistic analognes of ^/s/{2c^), v is the analogne of c, and n 
corresponds to the nnmber of extra dimensions. 


When b is not too small, the qnantnm-mechanical scattering phase of the wave associated 
with angnlar momentnm ih = mvb is 

1 


5 .= 


vh 


drV (r). 


(34) 


The intergral in eq. (^) converges only for n > 0, leading to 6b = —{bc/b)'^, where be — 
[Gn'mM/{hv)Y^'^ exactly corresponds to the dehnition in eq. (pISl). Therefore, the forward 
amplitnde (which is proportional to / db b6b) is hnite only for n > 2. For n = 0 one hnds 
the Conlomb singnlarity, and for n = 2 a logarithmic divergence. This is in agreement with 
eq. (^) and with the relativistic resnlts presented in sect. B 


A classical description of the process is valid if the qnantnm-mechanical nncertainties in 
the impact parameter b and in the scattering angle 6 are small with respect to their classical 
valnes. Using the Heisenberg principle, the qnantnm nncertainties are estimated to be 




(36) 


mv mvAb ’ 

where q is the momentnm in the direction orthogonal to the initial velocity v. The classical 
scattering angle is eqnal to the force {dV/dr) at the minimnm distance (r = b), times the 
collision time (b/v), divided by the momentnm (mv), 


b dVjb) 
mv"^ db 


(36) 
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The condition for the validity of the classical approximation {A6 -C 0, Ab -C b) implies 
\dV{b)/db\ > hv/b‘^. For a Coulomb-like potential V{r) = a/r, this condition is satished in 
the non-perturbative region a > fiv. In our case, it implies (for n > 0) 6 < 6c, with be ~ 
[GcmM/. The length scale b^ is a quantum scale {be ^ oo as h ^ 0) characteristic 
of the extra-dimensional scattering (it is undehned at n = 0). Therefore, at 6 > be, quantum 
mechanical effects cannot be neglected. 


In the classical process (for n > — 1), a projectile with an arbitrarily large impact param¬ 
eter b is deflected by a non-vanishing (although very small) angle, and therefore the cross 
section diverges at g = 0 (6 —cx)). In our case, for b > be, quantum mechanics sets in, and 
the quantum fluctuations of the scattering angle become larger than its classical value, ren¬ 
dering ill-dehned the notion of a trajectory. This is why the quantum scattering amplitude 
can remain hnite at g —0. 


It is also interesting to compare the result of eq. (h^ with the fully D-dimensional case, 
in which both gravity and the colliding particles propagate in the extra dimensions IQ. In 
this case eq. (|T^) becomes 


.4,1 = -2is / _ l 


(37) 


For g = 0, the integral is infrared dominated by large values of b and it is quadratically 
divergent (for any n). Therefore we hnd Aeikid —> 0) ~ b'^/q^, and we encounter the Coulomb 
singularity characteristic of long-range forces. On the other hand, in the case of particles 
living on the 3-brane, we have found (for n > 2) that the forward scattering amplitude is 
hnite and that, although the force does not have a hnite range, the integral in eq. m at 
g = 0 is dominated by distances of order be- Indeed, by working with localized particles we 
lose momentum conservation in the transverse direction, and moreover the Born amplitude 
becomes dominated by contact terms of ultraviolet origin. The fact that eq. (pTl) is not 
infrared dominated even for g —0 is a remnant of this property of the Born amplitude. 


Let us now come back to the relativistic elastic scattering under consideration. As we have 


seen in sect. |2.1| , there exist two different kinematic regions in which the eikonal amplitude 
has distinct behavior. 
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Region (i)\ y/s q 

This is the region in which the integral in eq. is dominated by the stationary-phase 
value of the impact parameter b = bg, where 


bg — be 




2(n+l) 


n+2 

2(n+l) 


Mi. 


(38) 


In this region b < be, and therefore the eikonal phase y = (be/b)"' is large (in units of h): we 
are in the classical domain. Indeed, the classical scattering angle is given by the derivative 
of the scattering phase with respect to the angular momentum L = y/sb/2, 


dx 2nV{n/2) Goy/s 
dL ~ 7r’"/2 


(39) 


In the limit n —>• 0, we recover the Einstein angle 6c\ = AGo^fs/b while, for n > 0, eq. (|H^ ) 
gives its higher-dimensional generalization^]. Notice that the relation between scattering 
angle and impact parameter in eq. (^) is equivalent io b = bg. 


Region (ii): q <^b~^. 

In this region the integral in eq. (|20|) is dominated by b of the order of (or slightly smaller 
than) be- This means that the eikonal phase y = {be/b)^ is of order one (in units of h) and the 
quantum nature of the scattering particles is important (although the exchanged graviton is 
treated as a classical field, and so quantum gravity effects are negligible), as we could have 
expected from the discussion of the non-relativistic analogue presented above. There is no 
one-to-one relation between impact parameter and scattering angle in region (ii), since the 
classical concept of a trajectory is not applicable. Moreover, notice that the relevant x never 
becomes much smaller than 1, and therefore we never enter the perturbative regime in which 
a loop expansion for the amplitude applies. Even though the interaction vanishes at 6 —> exo 
(where y —> 0), we never reach the Born limit. Even for q = 0, the scattering is dominated 
by 6 = 6c and not by 6 = oo, as in the Coulomb case. As explained after eq. (|37D , this 
result simply follows from the different dimensionalities of the spaces relevant for scattering 

^The Einstein angle defines the deflection of a photon by the static gravitational field of a mass-M 
particle at rest: Oe = iGiyiM/b. It is easy to see that by boosting 0^; to the center-of-mass frame, in the 
limit y/s ^ M, we obtain our expression of 0ci- 
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particles (living on a 3-brane) and exchanged graviton (propagating in the bulk), and it does 
not hold in the case of scattering of bulk particles. 


2.4 Corrections to the Eikonal Approximation 


The conditions required by our approximations are expressed in eq. (|). In order to assess 
the theoretical uncertainty inh erent to our calculation, in this section we estimate the size 
of the effects that we have neglected. We hrst discuss classical effects, for which general 
relativity suffices. We also include the quantum-mechanical effects of region (ii), but we 
leave quantum-gravity contributions for the next section. The neglected classical terms 
come from the approximations of small angle {—t/s -C 1) and of weak gravitational held 
{Rs/b <C 1). Associated with the subleading classical effects there is also the emission of 
gravitational radiation, leading to missing energy signals at high-energy colliders. 


Let us estimate these effects. Notice that in the classical region (i) where a/s ^ q ^ b 


-1 
C 5 


the two requirements coincide since h = hs, see eq. (|38|) . Here, scattering can be described in 
terms of classical trajectories. The equation b = bg represents a relation between transferred 
momentum and impact parameter 

2n+2 

(40) 


t 




D-- 


y2n+2 


b 


where we have neglected factors of order unity. This result represents an approximation 


in which one considers only a linear superposition of the gravitational shockwave helds 


generated by the two colliding particles []I^, thereby neglecting non-linear effects of their 
mutual interactions. In this computation, each held represents an exact solution to Einstein’s 
equations in the absence of the other held. Inclusion of the non-linear ehects will induce 
0{G^) corrections to eq. (^0]) , where p is an unknown power. By dimensional analysis the 
relative size of these corrections must be G^j^s^R/bPb<-^^'> r'U {Rs/by^'^^^\ Moreover since we 
expect that the evaluation of these ehects is perfectly perturbative, analyticity in s should 
hold, and therefore only even integers of the exponent p will appear. We conclude that 
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(41) 


eq. (0) will be corrected by a factor of the form 


1 + 0 




1 + 0 



This will induce a correction of the same size to the cross section, in the classical region 
(i). Notice that eq. (^H) corresponds to a 0{6^) correction term to the relation between 
scattering angle and impact parameter, eq. (|39|). In the scattering by a static potential we 
would get corrections already at 0{6‘^). Our different result follows from the relativistic 
nature of the process. 


In region (ii) where q the leading correction to the cross section will then have 

the form 

da.iw r. t ( 42 ) 


da 

dt 


eik 


dt 


t 


1 + - Zitbi) 


where Z is expected to be a constant of order unity in the classical region S> 1. In the 
region -C 1, the quantum nature of the scattering particles becomes relevant, and Z will 
have non-trivial behavior. However, admitting only hnite corrections at t 0 implies that 
the size of the relative correction to da/dt is at most 1/sb'l {Rs/b,r^+^. 


In summary, we have inferred that the leading classical gravity corrections to the eikonal 
amplitude are 



(43) 


and therefore suppressed and under control for the conditions we are applying, as expressed 
by eq. (1). 


The discussion outlined above is consistent with the analysis of ref. 0 where classical 
corrections to the eikonal are identihed by resumming an improved ladder series including 
a class of two-loop graphs, the so-called H-diagrams. After improving the Born term by the 
H-diagram, the eikonal phase is modihed to 


X 



1 + 0 



(44) 
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It is easy to check that altering eq. ( 0 ) according to the above equation produces corrections 
to the leading amplitude consistent with eq. (|^). Notice that the corrections to the eikonal 
become large for impact parameters comparable to the Schwarzschild radius, for which the 
production of black holes presumably sets inQ 

Gravitational radiation is associated with subleading classical effects. Emission of gravi¬ 
tons is signaled by the presence of an imaginary part in the H-diagram contribution to the 
eikonal phase |0] 

Im(x^) ~ , (45) 

br = . 

This absorptive term corresponds to a depletion of purely elastic scattering due to the emis¬ 
sion of gravitational radiation. When b br the probability of elastic scattering becomes 
exponentially small. The calculation of ref. 0 shows that the emitted radiation has typi¬ 
cal transverse momentum ~ 1/6. However the longitudinal momentum is distributed up to 
values of order ^/s. Therefore we expect that for b < br a signihcant fraction of the energy 
is radiated in the form of forward gravitational radiation. Notice that 6 < 6,. is typically 
outside our chosen kinematical regime, so that we expect only a fraction (6^/6)3^+2 ^ X 
of initial energy to be lost to invisible gravitational radiation. In the semiclassical region, we 
hnd {br/b)^'^~^‘^ ~ (—and therefore gravitational radiation 
for small-angle scatterings is not large. 

One may worry that important emission of radiation, already at 6 S> niay drastically 
reduce the naive geometric estimate of the black-hole production cross-section. Indeed if 
one simply interpreted eq. (^), a la Block-Nordsiek, as the number N of emitted gravitons, 
one would conclude that hard bremsstrahlung depletes all the energy, as soon as 6 is smaller 
than br- However this naive interpretation is manifestly inconsistent as the total emitted 
energy would quickly exceed ^/s. Indeed it is likely that the exponentiation of the H-diagram 

^The hypothesis that classical physics determines the black-hole production cross section has been criti¬ 
cised in ref. ||2^. However, we believe that the appearence of corrections of the form shown in eq. (Q) is a 
further indication that non-trivial classical dynamics emerges at 5 ~ Rs- 
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imaginary part is inconsistent for hard radiation. On the other hand for soft radiation with 
E 1/6, the Block-Nordsiek interpretation is likely to be reasonable. Thus if the number of 
soft gravitons is A^soft ~ (6r/6)^"’'"^ then the total energy lost to soft radiation at large angle 

is 



which becomes important only for b ~ Rs- 


2.5 Quantum-Gravity Effects 


On general grounds we expect that quantum-gravity contributions to the eikonal phase must 
be of the form Sx/x ~ with fc > 0 and \p given in eq. (^. If we work in the true 

transplanckian region {^/s S> Mp), these effects are subleading]^ since, as we have discussed 
in sect. 1, the classical scale Rs is much larger than the quantum scale Ap. 


However, as it will be clear in the following, collisions at the LHC can only barely sat¬ 
isfy the transplanckianity condition. Therefore, quantum-gravity corrections are potentially 
large, and they are likely to be the limiting factor of our approximations. This happens be¬ 
cause, although quantum-gravity effects are characterized by a length scale Ap smaller than 
the Schwarzschild radius Rs, classical corrections appear with a large exponent, see eqs. (¥T) 
and (1^), remnant of the geometric nature of the interaction, while this is not expected to 
be the case for the quantum-gravity corrections (9(Ap/6)^. 


However in order to be able to say more we need a quantum-gravity model or at least a 


framework. String theory is the only possibility at hand, and the transplanckian regime has 
been discussed in this context by several authors, see e.g. refs. 0. [n], H. Indeed, we will 


show that string corrections appear to be large for energies relevant to LHC experiments. 


"^It is interesting to compare the case of gravity with the case of vector-boson exchange which has been 
mentioned in sect. 1. For a vector, the eikonal phase is xv oc e^/6" = iRMgbR (we take h = c = 1). 
Notice that, as opposed to gravity, there is no factor of s enhancement |^. It is easy to realize that 
quantum corrections will also go like a power of e^/6". This is for instance the case for vacuum-polarization 
improvement of the vector propagator. Therefore, there is no impact parameter choice for which both 
Xv ^ 1 and quantum corrections are small, where the eikonal approximation can be useful. 
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In ref. 0] string corrections to eikonalized graviton scattering for type II superstrings were 
studied. Since we do not stick to a particular string realization of our brane-world scenario, 
we will limit ourselves to a qualitative discussion. Following ref. p9 we may imagine a 


realization of the brane-world in type I, where m = Q — n spatial dimensions out of the 10 
space-time dimensions are compactihed on a torus T"* with radius = 1/d^s, and 

Ms is the string scale. Of the remaining dimensions, n have a “large” compactihcation radius 
and 3 are the usual non-compact ones. The brane-world is realized by a (3 -|- m)-dimensional 
D-brane, in which m dimensions span the small manifold T™. In terms of Ms and of the 
string coupling gs, the 4 -|- n dimensional Planck scale M^ and the gauge coupling g on the 
brane are given by 


= 


TT 


9^ 


2 % ’ 


9s = 


(47) 


where 0 is the dilaton vacuum expectation value. Of course these relations will receive 
corrections of order unity in realistic models with broken supersymmetry and should be 


viewed, here and in the following, only as indicative estimates. Equations (^71) imply 


Me 


M, 


D 


2+n 


r 

47r 


(48) 


If we identify g with the gauge couplings of the SM, say for instance g = g^eak, we conclude 
that there should be a mild hierarchy between Mj:, and Ms- For instance, we hnd M^/Ms 
equals 1.7 and 2.8 for n = 6 and n = 2, respectively. A priori one could also conceive 
string models where M^, and Ms coincide, with a self-dual dilaton vacuum expectation 
value 0 ~ 0, corresponding to a strongly coupled string. In this case the weakness of the SM 
gauge couplings would have to be explained by some additional mechanism. 


As we will discuss in sect. |3)2, in order to observe transplanckian scatterings at the 
LHC, we have to consider values of M^, smaller than a few TeV. The existence of string 
theory at such a low scale is already limited by precision electroweak data and by the non¬ 
observation of new contact interactions at LEP. Exchanges of gravitons and massive string 
states give rise to effective operators involving SM helds. The leading effects come from 4- 
fermion operators of dimension 6 which, if present at tree level, give an approximate bound 
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Mg > 3 TeV 1^. This bound is in conflict with the working assumptions necessary to 
have observable transplanckian signals at the LHC, and therefore we have to assume a mild 
suppression of the dimension-6 operators. This assumption is probably not unreasonable, 
as there exist examples [^, where these effects vanish at tree level. In this respect 
the presence of supersymmetry at some stage in the construction may help suppress these 
contributions. Loop effects are within current experimental bounds even for Mu of order a 
few TeV, if the underlying quantum-gravity theory does not become strongly-interacting or, 
in a more operative sense, if the divergent loops are cut off at a scale slightly smaller than 
Md [^|. This allows for the weakly-interacting string scenario in which there is a small 


hierarchy between Mg and Mb, but it essentially rules out the case of strongly-coupled 
strings with Mg of a few TeV, where we would need the rather strong, and implausible, 
assumption that dimension-6 operators are suppressed to all orders in perturbation theory. 


With these limitations in mind we now consider string corrections [0. The basic result 
is that the eikonal elastic held-theory amplitude gets promoted to a unitary matrix acting 
non-trivially over a subspace of the string Fock space 


exp 




K = exp 





K, 

aX- 


(49) 


Here iL is a hermitian operator involving the creation and annihilation operators for the 
string oscillator modes. The above equation shows that, for b < bo, sl diffractive production 
of excited string modes takes place without any suppression. Basically what happens is 
that the two colliding partons are excited into string modes through multiple soft graviton 
exchange. 


The average mass of the string modes is roughly given by Mglbo/b)"'^^. This is smaller 
than the available energy y/s as long as {Rg/b)'^^^{\/^/b) < 1. This condition is clearly 
satished within our chosen kinematical regime, and therefore the produced particles are 
relativistic. The opening of an (exponentially) large number of new channels drastically 
depletes the elastic one. For instance in the example discussed in ref. 0, where gravitons 
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are scattered in type II superstring, the elastic amplitude correction factor becomes 


{K) = r2^+2(i _ iA)r2[l - i{n + 1)A], 


( 50 ) 


where A = (n/2)(6£)/6)"'+^ and where by {K) we indicate the expectation value of the 
operator K on the initial state. 


For b ':$> bo, the correction factor in eq. (^) becomes 


{K) ~ 1 + i 2 7 n(n + 1) 



(51) 


where 7 = 0.577 is the Euler number. Therefore, the effects of string excitations set in 
much before the string scale. Already at 6 ~ be, the terms modifying the gravitational 
eikonal phase implied by eq. are as large as n{n + As we have 

previously anticipated, these effects can be much more sizable than the classical corrections 
in eq. (|43|) , even if the string length is smaller than the classical Schwarzschild radius. 


For b -C bi), we hnd 


A| ~ (n + 



—7rn{n + 1 ) 



(52) 


showing that the absolute value of the correction factor decreases exponentially at b < b^, 
because of the appearence of the string excitations. 

This string effect is fairly general and we expect it in realistic models not to differ too 
much from eq. (|50|). However the sharp drop in the elastic cross section can be compensated 
by looking at the inclusive diffractive cross section. Since bjy < be we can use the saddle 
point approximation to evaluate the amplitude in transverse momentum space 

bp \ 
bs{q)) 

where ^eik is the eikonal amplitude in the stationary phase regime. Due to unitarity of the 
second factor, the total cross-section obtained by summing on the diffractive channels is just 


A{q) = Aeikiq) exp 


H 


(53) 
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the same as in the absence of strings, i.e. there is no exponential suppression. Notice that 
this inclusive cross section is consistently dehned in the limit in which the string modes are 
relativistic: at a hxed q all modes come out at essentially the same angle 6. The question of 
how well one can measure experimentally the inclusive cross section depends very much on 
the model, and in particular on the decay properties of the produced string modes. These 
states are excited by gravitons, so that they have the same color and electric charge of the 
incoming partons. If they are stable, the hnal state of the collision is given by two (massive) 
jet events. If they are unstable and their dominant decay mode involves bulk states, like the 
graviton, then the presence of missing energy would make it impossible to deduce the original 
center-of-mass energy and to make sure that only events at a given center-of-mass energy 
are selected. On the other hand if they dominantly decay on the brane this reconstruction 
might be possible. Notice that the coupling to bulk modes (closed strings) and to brane 
modes (open strings) are proportional to and gs, respectively. Therefore the decay to 
brane modes, when possible, dominates over bulk decays in weakly-coupled string models. 


We also want to comment on the string effects arising in the regime Mo < ^/s < 
Ms/gl ||2^ . In this regime, although the energy is transplanckian, Xs > Rs holds. Then 
the large angle, small impact parameter effects are dominated by string physics instead of 
classical gravitational effects. In this regime one should use the string version of the eikonal 
phase 0. The main effect is that Imy ~ g^sa' exp{—b‘^/2a' In sa'), corresponding to inelas¬ 
tic production of very excited strings (string balls) dominating the scattering for b < Xs- 
The corresponding cross section is ~ a' for Ms/gs < \fs < Ms/gs- At y/s Ms/gh we 
have that a' ~ consistent with the string balls becoming indinguishable from ordinary 
black holes . The onset of inelastic string production marks the end of the validity of our 
description, since any reminiscence of classical gravitational dynamics is lost. 


Notice that in the context of strings there exists another approach to calculate hxed- 
angle scattering in the regime y/s ^ Ms [^. This different approach leads to an elastic 
amplitude that vanishes at small hxed angle like exp{—9y/a's) |^. Recently this approach 
was applied m to the TeV string scenario. The conclusions of ref. il would seem to 
disagree with those of ref. [0, based on the string eikonal, according to which the scattering 
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at y/s S> Ms is well described by (semi)classical general relativity. However, as explained in 
ref. [^, the Borel resummation leading to the result exp{—9^/a's) applies without relevant 
corrections only in a limited range of energies y/s < Ms{lia . Indeed for energies such 
that As ^ Rs the two approaches give results that are fairly consistent with each other K 


So we believe that in the transplanckian regime our simple approach is valid. Notice, that 
also in the approach of ref. 0 the elastic amplitude has an exponential suppression, different 
than that in refs. [^, |^, and due to gravitational radiation and diffractive string production. 
Moreover the discussion of ref. 0 shows that, even if the elastic channel is suppressed, the 
total cross section at hnite angle is large and grows with y/s. Notice also that if y/s is large 
enough the average mass of the diffractive states eventually exceeds both Ms and M^, so it 
is possible that the diffractive channel becomes at extremely high energies the production of 
a pair of small black holes. 

Next, we would like to comment on the possible effects of the modes corresponding to the 


brane motion, the branons [^. The hrst comment is that these modes do not necessarily 
exist. If the SM helds are localized at some hxed point, as in orbifold constructions, then it 
is consistent to assume that the corresponding “brane” is not a dynamical object, and could 
also have vanishing tension. So when the SM lives at a hxed point there is no effect (apart 
from a “trivial” reduction of the phase space of bulk gravitons). 

Two types of effects can be associated with brane dynamics. First, the presence of a 
brane tension r = /i^ gives rise to a gravitational held, which at a tranverse distance y from 
the brane vanishes like {Rr/y)^~‘^, where the gravitational radius Rr is 

r(n/2)2^+Vt 


R^-^ = 


T 


(54) 


^(n — l)(n + 2) M'^^ 

The wave function of a bulk graviton with transverse momentum kx is expected to receive 
a correction of order {RT-kT)"'~^ from the presence of this background. For a given impact 
parameter b, the tree-level scattering amplitude is precisely saturated by the exchange of 
graviton KK modes with mass ~ 1/6. Therefore we expect the background metric to correct 
the eikonal phase by a relative amount ~ (i?^/6)"'“^. For n > 2, at large enough impact 
parameter we can neglect the brane radius as much as we can neglect the Schwarzschild radius 
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Figure 3: Ladder and cross-ladder Feynman diagrams contributing to elastic scattering of <F 
particles. The dashed lines represent the exchange of brane excitations. 


Rs- This gravitational brane effect becomes less important as the tension is decreased. 


The second class of effects is associated with the exchange of branons. Branon exchange 
becomes more important the smaller the tension. Notice hrst of all that at the classical level 
branons are not excited by the gravitational shockwave of a fast moving particle, i.e. the 
induced metric on the brane remains of the Aichelburg-Sexl kind. At the quantum level, 
branons can be exchanged in pairs between the colliding particles. The exchange of a pair 
of branons is very much equivalent to the exchange of a “composite” graviton. In order to 
get an estimate of the size of the effect, we may generalize the graviton eikonal including the 
ladder and cross-ladder diagrams with two-branon blobs (see £g. The eikonal phase gets 
modified into 

= (v) 

In the case of n = 6 the effect is equivalent to a renormalization of Gn, but in all other cases 
it can distort our result from pure graviton exchange. 


Including both effects from brane excitations, we can write the relative correction to the 
elastic amplitude (in the semiclassical region) as 

n —2 Q 

where we have dropped all numerical factors. This equation shows that for the interesting 
case 2 < n < 6 and for the natural choice /i ~ Mp, both branon effects are small in the region 
of small angle and large s. Notice that for a 3-dimensional D-brane arising by wrapping a 
(9 — n)-brane over the 6 — n dimensions of radius a/^ along the lines of the example above. 


b — n 
n+1 


(56) 
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the tension is given by 


r = 


(27r) 2Q;'^ 


(57) 


Then at weak conpling we have Rr < v a' and the effects of the first class are very small. 
Effects of the second class are also small as long as b > and n < 6. For n = 6, 
the contribntion from branon exchange leads to the pecnliar resnlt Xbranon = Xgrav/5, see 

eq. ®)Q 


To snmmarize, qnantnm-gravity corrections are potentially very sizable because, although 
they involve smaller length scales, they do not appear with large exponents, as in the case 
of the semiclassical corrections. Of course, they cannot be computed in the absence of a 
complete quantum-gravity theory. We have estimated their effects in the context of string 
theory, borrowing the results of ref. (which are however strictly valid only for graviton- 
graviton scattering). For b < bjj, a diffractive production of string states occurs, depleting the 
elastic scattering cross section, although it may still be possible to recover the gravitational 
scattering properties by studying inclusive cross sections. This gravitational excitation of 
string modes looks like a very promising way to investigate, in a rather clean channel, string 
effects at high-energy colliders. As the impact parameter is further decreased and becomes of 
the order of \s = \/^, we are entering the regime of head-on collisions between string modes, 
with inelastic production of Regge excitations eventually leading to multi-string states Hi- 


In the following, we will focus on gravitational transplanckian collisions, neglecting 
quantum-gravity effects. One should keep in mind that such (at the moment incalculable) 
effects may give significant modifications to the signal discussed below. Although potentially 
disrupting to the predictivity of the elastic channel, these new-physics effects are of course 
very interesting from the experimental point of view. 


3 Phenomenology in the Transplanckian Regime 


^We thank A. Strumia for stressing this result of the physics of D-branes. 






3.1 Signals at the LHC 


At the LHC, the observable of interest is jet-jet production at small angle (close to beam) 
with large center-of-mass collision energy. The amplitudes derived in previous sections are 
applicable for the scattering of any two partons. The total jet-jet cross-section is then 
obtained by summing over all possible permutations of initial state quarks and gluons, using 
the appropriate parton distribution weights and enforcing kinematic cuts applicable for the 
eikonal approximation. 


Dehning s and t as Mandelstam variables of the parton-parton collision, we are interested 
in events that have ^/I/Md S> 1 and —i/s -C 1. We can extract from the jet-jet invariant 
mass Mjj = \/S, and i from the rapidity separation of the two jets —i/s = l/(l-|-e^'^), where 
Arj = Tji — r] 2 - The variable Ar/ is especially useful since it is invariant under boosts along 
the beam direction, and it is simply related to the 9 scattering angle in the center-of-mass 
frame: 


A ?7 = In 


1 -I- cos 6 


1 — cos 9_ 

Therefore, the kinematical region of interest is defined by the equivalent statements 


(58) 


^7] OO 0 —> 0 



(59) 


In computing the cross-sections we use the CTEQ5 ||3^ parton-distribution functions. 
We have also compared results with GRV [^] and find little difference. We evaluate the 
parton-distribution functions at the scale = 6“^, see eq. (|3^) , if g > h~^ and 
otherwise (g^ = —i) 


We require that both jets have |? 7 | < 5 and pt > 100 GeV from conservative detector 
requirements. The SM di-jet cross section is computed using Pythia |^, ignoring higher- 
order QGD corrections, and the gravitational signal is calculated using our own VEGAS 
Monte Garlo integrater. For simplicity we are dehning the background as the jet-jet cross- 
section from QGD with gravity couplings turned off, and the signal as the jet-jet cross-section 
from the eikonal gravity computation with QGD turned off. In reality, SM and gravity 
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contributions would be simultaneosly present. Nevertheless, there is no interference between 
the leading QCD contribution in the limit t/s —0 and the gravitational contribution, since 
graviton couplings are diagonal in color indices. However, terms 0{aw) are generated by the 
interference between gravitational and Z /7 contributions. At order the di-jet cross 

section receives contributions from the square of the amplitude obtained by exchanging 1 
gluon and any number of gravitons (resummed in the eikonal approximation) and from the 
interference between the gravitational amplitude and resummed diagrams with exchange of 
2 colored particles and any number of gravitons. In the small regions of parameter space 
where the QCD and gravity contributions are comparable, these terms, which will not be 
computed in this paper, should be taken into account. We will not make any detailed claims 
about that region here. 


Although we are working in the transplanckian regime, we are considering processes with 
low virtuality {Q^ ~ 1 / 6 ^ M^) and, as long as no new physics appear at momenta smaller 

than Q, the SM dynamics give a reliable estimate of the background. In the context of 
string theory, the QCD background has to be interpreted as the leading contribution from 
exchange of open strings, in the limit in which the transfer momentum is smaller than the 
mass of the Regge excitations. 


The jet-jet process is best described by the two kinematic variables Mjj and At]. Although 
experimental searches for new physics in the jet-jet observable should investigate the full 
range of these two kinematic variables, our goal is to study the region in which the eikonal 
computation is valid, where we can make reliable theoretical predictions. This leads us to 
the region of large At] and large Mjj. 

Let us begin by studying the parton differential cross section signal distribution as a 
function of the rapidity separation Arj for a hxed \/§ = Mjj, 


da 

dAf] 


Tiblse^^ 

{1 + e^^f 


Fn{y)\\ 


(60) 


where y 


This distribution has a series of peaks whose maxima and minima 
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are determined by the values of Arj that satisfy the equation 

y d\Fn{y)\ 


\Fn{y)\ dy 


bcVs 


( 61 ) 


y= 




These peaks arise from the oscillations of the function \Fn\ around its asymptotic value 
given in eq. (^), and characterize the transition region between the “classical” (small Arj) 
and “quantum” (large Arj) regimes discussed in sect. Therefore, we cannot solve eq. (|^) 
using the asymptotic expressions given in eq. (^^ or eq. (p5[) . Approximate solutions of 
eq. (^) can be found if the peaks are located at sufficiently large values of Arj, since in this 
case we can neglect with respect to 1 in the left-hand side of eq. (|6l|). Then, eq. (|61|) 
becomes a function only of the variable y. For our consideration, this approximation is 
adequate for determining at least the hrst peak. We hnd that the hrst peak of the Ap 
distribution, for a hxed value of the two-jet invariant mass Mjj, is given by 


Ar/(P"ak) = In 




n 


M. 


(62) 


D 


The coefficients can be obtained by numerically solving eq. (l6l|) and are given by ^ 2 , 3 , 4 , 5,6 = 
0.8, 0.9, 1.0, 1.2, 1.3. 

These peaks are a characteristic feature of the higher-dimensional gravitational force, 
and correspond to the diffraction pattern of the scattered particles. As discussed in sect. 
higher-dimensional gravity, albeit leading to a force with inhnite range, dehnes a length scale 
[be) in the quantum theory. A diffraction pattern emerges when q ^ b~^ corresponding to 
the interference of the scattered waves in the region b ~ be- In the Coulomb case (u = 0), 
such a scale does not exist and therefore no diffractive pattern is produced. 

The di-jet differential cross section dajj/d\A7]\ is plotted in £g. || for u = 6, Mjj > 9TeV 
and Mr, = 1.5 TeV and 3TeV. Since the parton-distribution functions decrease rapidly at 
higher Mjj, the value of Mjj in eq. (|6lD is well approximated by 9 TeV in this example. The 
hrst peaks are then calculated from eq. (|6^ to be at Ay = 5.5 {Mr = 1.5 TeV) and 3.7 


{Mr = 3TeV). The Monte Carlo integrated distributions have peaks that agree well with 
these numbers. 
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Figure 4: The di-jet differential cross section dajj/d\Ari\ from eikonal gravity for n = 6, 
Mjj > 9TeV, when both jets have \ri\ < 5 and pt > 100 GeV, and for = 1.5 TeV and 
3 TeV. The dashed line is the expected rate from QCD. 

The same differential cross section is shown in £g. ^ for different values of n and for 
Md = 1.5 TeV. For n = 2 the first peak is partly hidden by the logarithmic divergence for 
t —>■ 0, and no structure after the hrst peak is visible. As n increases, the peaks become 
more evident. The study of the peak structure could be a feasible experimental technique 
to measure the number of extra spatial dimensions n. However, such a study can only 
be pursued with a full detector simulation, taking into account the rapidity and jet mass 
resolutions. 

Since the two jets are experimentally indistinguishable, we have used |A? 7 |, instead of Ar/, 
as the appropriate kinematical variable to plot. This means that the experimental signal 
considered here contains also contributions from scattering with large and negative Ar/, 
which corresponds to partons colliding with large momentum transfer and retracing their 
path backwards. For the background, these effects are calculable and taken into account. 
However, the theoretical estimate of the signal at negative At] lies outside the range of 
validity of the eikonal approximation. Nevertheless, this is expected to be negligible and can 
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Figure 5: The di-jet differential cross section dajj/d\Ari\ from eikonal gravity for M^, = 
1.5TeV and n = 6 (solid line), n = 4 (dotted line), and n = 2 (dashed line). We require 
Mjj > 9 TeV and that both jets have \ri\ < 5 and px > 100 GeV. 


be safely ignored. Indeed, by dimensional arguments, the differential cross section for t —>■ —s 
is estimated to be da/dt ~ nRg/s, see eq. (^Ol). When compared to the signal at small t 

in+2p 

{da/dt ~ vr6^), this gives a correction 0[Rg/{sb"/)] ~ C>[(M|,/s) "("+!>], which is smaller than 
the terms we have neglected in our calculation, see eq. (^31). Notice that this complication is 
unavoidable in the presence of identical colliding particles, since the amplitude is symmetric 
in t and u {u = —s — t). 


Recalling our discussion of the corrections to the eikonal amplitude given in sect. |2]J, we 
notice that semiclassical corrections are indeed small. In the kinematical region of interest. 


the second term in eq. (^) amounts to less than 1% for Mjj = QMb and n < Q and less 
than 5% for Mjj = SM^ and n < 6, while the first term amounts to about 5% at Ap = 3. 
Gravitational radiation at Ap = 3 and Mjj = 3 Md gives corrections of about 7% for n = 2 
and 5% for n = 6, see eq. (^). Quantum gravity can however severely affect our signal. If 
the corresponding corrections had the form {\p/hcY, then they would amount to 5% (6%) 
for Mjj = QMp and n = 2 (n = 6) or to 20% (9%) for Mjj = SMp and n = 2 {n = 6), 
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and they would be under control. However, it is rather likely that the dynamics taming the 
ultraviolet behavior of gravity sets in at energies lower than Md, enhancing these corrections. 
This is indeed the case of weakly-coupled strings, where in the previous estimate Ap should 
be replaced by A 5 , giving an enhancement of For instance, if the diffractive string 

production discussed in sect. takes place, it will deplete the elastic channel at impact 
parameters smaller than fcp, which is typically very close to be, for values of Mp relevant 
for the LHC and for weakly-coupled strings. In this case, string production could be the 
discovery process at the LHC. 


As is apparent from fig. the background increases with Ap faster than the signal, and so 
the best signal to background ratio is found at the smallest Ap. However, small Ap are not in 
accord with the eikonal approximation. To stay within the acceptable kinematic range of the 
eikonal approximation we impose the condition Ap > 3, which is equivalent to —t/s < 1/21 
and 9 < 25°. We then select a maximum value of Ap in order to reject the background. In 
our study we make a simple universal choice Ap < 4, but in practice one could choose the 
maximum Ap using only the criterion of retaining enough events to be detected. We want to 
stress that this universal range of Ap is chosen just for illustration, since its optimal choice 
depends on Mp. In particular, if Mp is below about 2 TeV, visible signals certainly can be 
obtained with a more stringent lower limit on Ap, therefore investigating a region of smaller 
scattering angles, where the eikonal approximation is even more accurate. 


Furthermore, with the current computation technology of QCD, the dijet rate in our 
kinematic conhguration is not a precision observable to compare with theory. This is espe¬ 
cially true at large rapidity separation. However, our gravity signal can swamp the QCD 
expectation most notably at smaller rapidity separation, where the hard-scattering QCD 
rate is smaller and more reliably known. That is partly the reason why we choose for some 
plots to illustrate the gravity signal in the modest interval 3 < |Ap| < 4. 


At very large rapidity separation, the QCD dijet observable may be dominated by BFKL 
dynamics |^, in which case a decorrelation of the azimuthal angles of the two jets ensues . 
Such a decorrelation is not expected in the color-singlet gravity signal. It may even be 
possible to use this to distinguish QCD from gravity. Similarly, the gravity signal is not 
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Figure 6: Total integrated di-jet cross-section for 3 < \Ari\ < 4, n = 6, and Mjj > 
when both jets have \ri\ < 5 and px > 100 GeV. Lines are plotted for M^, = 1.5 and 3 TeV. 
The eikonal approximation is reliable only where Mjj/Mx, S> 1. The expected QCD rate is 
given by the dashed line. 

expected to hll the central region with soft gluons, and the lack of this central jet activity 
could be an additional discriminating tool. 

After having selected the range of 3 < |A? 7 | < 4, we now show in £g. ^ the cross-section 
as a function of minimum jet-jet invariant mass cut for Mx, = 1.5 TeV and 3 TeV. We plot 
results for all Mjj > Mx,, but we recall that the eikonal approximation is valid only for 
Mjj/Mx) S> 1. This plot shows the important feature that the signal cross-section is flatter 
in Mjj than the background. This enables better signal to background for larger Mjj cuts, 
which is the preferred direction to go for eikonal approximation validity. Therefore, one 
should make the largest possible Mjj cut that still has a countable signal rate for a given 
luminosity. 

Finally, in £g. ^ we plot the total integrated cross-section as a function of Md for Mjj > 
SMjx (left panel) and Mjj > 6Mx (right panel). We also have required 3 < IAt^I < 4. The 
two solid lines correspond to u = 6 (upper line) and n = 2 (lower line). They are not far 
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Figure 7: Total integrated di-jet cross-section for 3 < IAt^I < 4, when both jets have \ri\ < 5 
and pt > 100 GeV, for (left panel) Mjj > SMd and (right panel) Mjj > 6M/). In each plot, 
the upper solid line is for n = 6 and the lower solid line is for n = 2. The dashed line is the 
expected QCD rate. 


separated in this log plot. Figure ^ demonstrates several universal features. First, larger 
Mjj/Mo corresponds to a larger ratio of gravity signal to expected QCD rate, and to a more 
reliable applicability of the eikonal approximation. Larger Mjj/Mj:, also means less range of 
Mjj probed at the LHC since parton luminosity drops rapidly as Mjj approaches the 14 TeV 
collider limit. In this paper, we will restrict our considerations to Mjj/Mjj > 3 and so elastic 
scattering with Md as large as 3.5 TeV can be tested with more than 10 events in 100 fb“^. 
If we adopted the more conservative constraint Mjj/Mr, > 6, the range of Mr, where large 
and calculable signals are expected reduces to 1.8 TeV. 

Before concluding this section on the LHC signatures, we recall that the eikonal jet- 
jet observable analyzed here is not the only transplanckian non-SM process. Black-hole 
production is also expected to be very large with a parton-parton production cross-section 
given approximately by arn — 7ri?|. As discussed earlier, this is just a dimensional-analysis 
estimate and reflects only the expectation that all collisions with impact parameter b < Rs 
get absorbed into a black hole. 

The description of the scattering process in terms of black-hole production is expected 
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Figure 8: Black-hole production cross-section as a function of (left panel) minimum black- 
hole mass and (right panel) the gravity scale M^,. 


to be reliable only for black hole masses Mbh well above the D-dimensional gravity scale 
Mb). In the left panel of £g. ^we plot the production cross-section as a function of minimum 
black-hole mass for M^, = 1.5 TeV and = 3 TeV. Although the lines are extended down 
to = Mo, the reliable region is only Mbh/Md S> 1. The signatures of black holes 

are spectacular high-multiplicity events with almost no SM background 0. Only a few 
events at very high invariant mass are needed to identify a signature. The right panel of 
fig. I shows the cross-section for black-hole production for Mbh > SMn and Mbh > 

With an integrated luminosity of 100 fb“^, a minimum of several events can be achieved for 
Mjj ~ 3.5 TeV, if we require Mbh > SM^j, or for Md ~ 1.8 TeV if Mbh > These 

values of Mj:, are nearly identical to what one can reach with the jet-jet observable discussed 
above. The combined search for black-hole production and small angle jet-jet events will be 
important for a full experimental characterization of the extra dimensions and for our ability 
to determine the underlying parameters of the theory. 


37 





3.2 Other Future Colliders 


In addition to Tevatron upgrades and LHC, there are many colliders envisioned for the 
future. These include variants on the e’''e“ linear collider theme, the muon collider, and the 
very large hadron collider (VLHC). Observables associated with gravity are very sensitive 
to collider energies. The search for phenomena of almost all beyond-the-SM schemes is 
interested in higher and higher energies, but gravity is especially dependent on this progress. 
The reason is that either the rate increases at a high power of the collision energy if s -C 
Ml, {e.g., a ~ for jet plus graviton signatures); or, well-chosen observables 

become calculable if s 3> Ml, {e.g., in the case of eikonal gravity jet-jet production). Such 
calculability is taken for granted in many other beyond-the-SM frameworks. Either way, the 
highest energies attainable are important for increased understanding if low-scale gravity is 
correct. 


We would therefore like to briefly discuss transplanckian collisions in the context of two 
of the highest energy colliders being considered in the not-too-distant future. The hrst is 
the VLHC [^, whose center-of-mass energy for proton-proton collisions is envisaged to be 
between 50 TeV and 200 TeV. Obviously the higher energy of this post-LHC machine will 
be much more probing. Much of the analyses that we have described for the LHC would be 
applicable for VLHC analyses. By going to higher center-of-mass energies, we can gain in the 
range of sensitivity to M^ and, more importantly, we can study kinematical regions where 
quantum-gravity effects are expected to be smaller and the eikonal approximation becomes 
more trustworthy. For instance, by analysing dijet hnal states at the VLHC with invariant 
masses larger than 150 TeV, we can make quantum-gravity corrections (proportional to 
A|/6^), for a hxed value of the string scale, about a factor of 200 (for n = 2) or 6 (for n = 6) 
smaller than their expected size in LHC experiments. 

The other collider discussed at the highest currently imagined energy is CLIC , which 
is a two-beam linear collider design. Energies as high as ^/s = 10 TeV are being 

considered. Although such high-energy machines will have some spread in luminosity as a 
function of the center-of-mass energy, one still expects that approximately 1/4 of the beam 
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luminosity will be within 1% of designed center-of-mass energy. With large gravity-induced 
cross sections (~pb or more) for all Mn/'/s -C 1, combined with relatively low background 
(~fb) and impressive luminosity design goals of 1 ab“^, CLIC could certainly improve, verify 
and even discover the gravitational origin of beyond the SM signatures. 


The maximum values of that can be studied at CLIC in the gravitational deflection 
process (Bhabha scattering) are uniquely determined by the condition for validity of the 
eikonal approximation. For comparison, the condition equivalent to the one chosen in our 
LHC analysis is < \/s/3. Therefore, a linear collider with ^/s = 10 TeV could probe a 
parameter-space region, where theoretical calculations are reliable, which is very similar to 
the one that can be studied by LHC. However, the cleaner e’''e“ environment offer several 
advantages for precision tests and parameter determinations. 

In £g. 1^ we give one example of the signal and background distributions in the scattering 
angle 6 for Md = 2 TeV and a/s = 10 TeV. Here 6 is the scattering angle of the electron in 
the Bhabha process, and 0 = 0 indicates the electron going down the beam pipe undeflected 
by the collision. Notice that the SM background is completely insignihcant as long as we 
exclude a small region around 0 = 0. This allows experimental studies of the cross section 
in a much more forward region than what is possible at the LHC. 


The signal angular distribution is characterized by the peak structure encountered before, 
with minima and maxima described by the equation 


2 V 2/ \Fn{,y)\ dy 

The hrst peak is approximatly given by 

^(peak) ^ 


(63) 


(64) 


where a„ is a numerical coefficient with value 02,3,4,5,6 = 0 - 9 ) 1-2, 1-1, 1.1, 1.0. The n = 6 
line (solid line) in fig. ^ is again the most telling line, since a careful measurement of the 
ups and downs of da/d9 would be hard to reproduce in another framework. Given the good 
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Figure 9: Angular distribution of the e'''e“ —> e’''e“ signal for = 2TeV and ^/s = 
lOTeV. The solid line is for n = 6, the dotted for n = 4 and the dashed for n = 2. The 
almost-vertical dash-dotted line at the far left of the hgure is the expected Bhabha scattering 
rate from the Standard Model. 

energy and angular resolutions that can be achieved at a linear collider, CLIC could perform 
a much more precise study of the peak structure than what is feasible at the LHC. 

Another advantage of CLIC is that unlike a pp collider, the scattering particles are 
distinguishable. This means that the calculable limit of t —0 (small angle scattering) is 
unambiguously identihable. Recall that in the pp collider case, the case of partons glancing 
off each other at t > 0 was not distinguishable from partons bouncing backwards at large 
momentum transfer i —s. We argued that these non-calculable large momentum transfer 
contributions are negligible, but CLIC can test that assertion. 


3.3 Comparison among Gravity Signals at the LHC 

Let us conclude by comparing the various signals and experimental strategies for the dis¬ 
covery and the study of higher-dimensional gravitational interactions at the LHC. The two 
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relevant parameters are the fundamental Planck mass Md and the number of extra spatial 
dimensions n. We are tacitly assuming that the new dynamics of quantum gravity does not 
introduce a new mass scale smaller than Md, or new phenomena in the relevant kinematical 
region. The searches for gravity signals at the LHC are classihed in terms of three different 
kinematical regions. 

Cisplanckian Region. This is the region in which the center-of-mass energy of the par- 
ton collision is smaller than the Planck mass, a/^ -C Md. The theory can be described 
by an effective held-theory Lagrangian with non-renormalizable interactions, and the seach 
for new contact interactions is an adequate experimental tool. However, the relations be¬ 
tween the coefficients of the contact interactions and the fundamental gravity parameters 
is model dependent. In this region, the search for jet plus missing energy events is partic¬ 
ularly interesting, because graviton emission can be reliably calculated with a perturbative 
expansion. The region |Q] in which we can trust the effective theory for graviton emission 
is Md > 3.8 TeV (for n = 2) or Md > 4.8 TeV (for n = 4). LHC is sensitive to the 
graviton signal up to Md <8.5 TeV (for n = 2 and an integrated luminosity £ = 100 fb“^) 
or Md < 5.8 TeV (for n = 4 and C = 100 fb“^). In this window, LHC can perform a 
quantitative test of high-dimensional gravity in the cisplanckian region. 

Planckian Region. This is the intermediate region \/^ ~ Md, where no experimental 
signals can be predicted without knowledge of the quantum-gravity theory. Signals from 
graviton emissions or elastic gravitational scattering are also present, but their rates cannot 
be reliably computed. New and unexpected phenomena will be the key to understand the 
underlying dynamics. Of course this is the region which will eventually yield the crucial 
experimental information. 

Transplanckian Region. This is the region discussed in this paper, characterized by 
a/s S> Md- The two experimental signals of interest in the transplanckian region are di-jet 
events (from elastic parton scattering) and black-hole production. Unfortunately, in the case 
of the LHC we are rather limited in energy for both processes. The actual conditions for 
calculability we have taken are > QMd (in which the transplanckian region extends up 
to about 1.8 TeV) or the more optimistic y/S > 3 Md (in which the transplanckian region 
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extends up to about 3.5 TeV). In both cases we are only marginally inside the transplanckian 
region since a complete separation between the quantum-gravity scale Xp and the classical 
Schwarzschild radius Rs has not been fully reached. Therefore quantum-gravity contri¬ 
butions can potentially modify in a signihcant way both elastic scattering and black-hole 
formation. In our analysis, we have assumed that such contributions are small. 

While the black-hole production cross section can only be estimated by dimensional 
analysis, the elastic cross section in the small-angle region can be computed as a perturbative 
expansion over controllable parameters. The elastic cross section is larger than the one 
for black-hole production. The observation of a cross section at hnite angle growing with 
a power of s would be a clean signal that the high-energy dynamics of gravity has been 
detected. Given the highly characteristic events from black-hole evaporation, we expect 
negligible backgrounds for the black-hole events. As we have shown in this paper, the QCD 
background for di-jet events can be overcome by studying the distributions in di-jet invariant 
mass and rapidity separation. Therefore, the di-jet signal from the gravitational deflection 
of partons is not limited by the background. Its characteristic distributions can be used 
to test the gravitational nature of the interaction and to determine the parameters of the 
underlying theory. 
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Appendix 


In this appendix we collect some mathematical formulae which are useful to reproduce com¬ 
putations presented in this paper. The Bessel functions Jn{x) have the following series 
expansions 


OO 


Jn{x) = Y. 

k=0 


kW(k + n + l)’ 


(65) 
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( 66 ) 



Some useful integrals are 


(2n + l)^ 


4n^ — 1 
8x 


sin 


X — (2n + 1)^ 



f (Ty e^^-yfiv) = f dy y'^ J’^_i{xy)f{y), 

J x^ JO 

for a generic function /, and with x = \x\, y = \y\. 

r dx =--T 

Jo ^ '' a \ nj 

/•e ^a+1+71 

/ dx x^ e*"'" =-[1 + 0 (e'^)] for a > -n - 1. 

Jo n 


(67) 
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